Four-body calculation of proton- 3 He scattering 
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The four-body equations of Alt, Grassberger and Sandhas are solved, for the first time, for proton- 
3 He scattering including the Coulomb interaction between the three protons using the method of 
screening and renormalization as it was done recently for proton-deuteron scattering. Various realis- 
tic two-nucleon potentials are used. Large Coulomb effects are seen on all observables. Comparison 
with data at different energies shows large deviations in the proton analyzing power but quite rea- 
sonable agreement in other observables. The effect of nucleon-nucleon magnetic moment interaction 
and correlations between p-d and p- 3 He analyzing powers are studied. 
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Modern calculations of light nuclear systems A < 12 
are essential to our understanding of the force models 
that have been developed to describe how nucleons in- 
teract at low energies 0, 0|- Of these nuclear systems, 
the four-nucleon (47V) system is particularly important 
because it gives rise, experimentally, to the simplest set 
of nuclear reactions that shows the complexity of heav- 
ier systems and the Coulomb interaction manifest itself 
in new ways relative to what is observed in the three- 
nucleon (3iV) system. Theoretically it is also impor- 
tant because with powerful numerical techniques and fast 
computers one can calculate not only bound state proper- 
ties Q but also scattering observables 0, d, d, 0, B S E3] 
for a number of elastic, transfer and breakup reactions 
that place new challenges to our understanding of the 
underlying force models. The importance of scattering 
calculations also has to do with the possibility to probe 
states in the continuum associated with specific reso- 
nances, states of higher angular momentum than corre- 
sponding bound states, effects that depend on the spin 
orientation of the projectile and/or target, and threshold 
effects on the observables, among others. 

While the three-nucleon system has been extensively 
studied [Til . Hit ] through neutron-deuteron (nd) and 
proton-deuteron (pd) elastic scattering and breakup ex- 
periments, exact calculations using realistic force models 
as well as interactions derived from Effective Field The- 
ory were restricted, for a long time, to the nd system 
due to limitations in including the Coulomb force in the 
description of pd scattering beyond low energy pd — > pd 
and pd <-> 7 3 He calculations Tj|, 14 1 in the framework of 
the variational hyperspherical approach. The situation 
has now changed due to the work of Refs. [TBI, where 
calculations of pd — ► pd, pd — ► ppn, pd <-* 7 3 He, 7 3 He — > 
ppn, e 3 He — ► e'pd, and e 3 He — * e'ppn were performed 
at energies ranging from 1 MeV in the center of mass 
(cm.) system to the pion production threshold. The 
work is based on the solution of the momentu m-sp ace 
Alt, Grassberger and Sandhas (AGS) equations [17| to- 
gether with the screening and renormalization approach 
lH, ISi 20] for the Coulomb interaction leading to the re- 



sults of observables that are independent of the screening 
radius, provided it is sufficiently large. 

In the present manus crip t for the first time we extend 
the method of Refs. fl5llld| to thep- 3 He elastic scattering 
using the four-body AGS equations 2l|. The aim is to 
bring the 47V scattering problem to the same level of un- 
derstanding in terms of the underlying two-nucleon (27V) 
forces as already exists for 37V, which means that cal- 
culations are carried out without approximations on the 
27V transition matrix (t-matrix) like in Ref. [f| or limita- 
tions on the choice of basis functions as in Refs. [22I l23j|. 
Therefore, after partial wave decomposition, the AGS 
equations are three-variable integral equations that are 
solved numerically without any approximations beyond 
the usual discretization of continuum variables on a fi- 
nite momentum mesh. The results we present here are 
converged vis-a-vis number of partial waves and momen- 
tum meshpoints as well as the value of the screening ra- 
dius of the Coulomb potential. These calculations are 
also an extension to p- 3 He of the work already developed 
for n- 3 H [24]], and were presented for the first time in 
Ref. [25j]. Our work follows the work of Refs. (JEHU, 
but with greater number of 27V, 37V, and 47V partial 
waves in order to get fully converged results for the spin 
observables and with various 27V potentials. The advan- 
tage of the present work is that it is easier to extended to 
inelastic reactions and to use with nonlocal interactions. 

Our description of 47V scattering is based on the sym- 
metrized four-body AGS equations given in Ref. 24 1 
where the solution technique is discussed in detail. In 
order to include the Coulomb interaction we follow the 
methodology of Refs. [la, [l6| and add to the nuclear pp 
potential the screened Coulomb one wr that, in configu- 
ration space, is given by 



!» H (r)^(r)e- W 



(1) 



where w(r) — a e /r is the true Coulomb potential, a e ~ 
1/137 is the fine structure constant, and n controls the 
smoothness of the screening; n — 4 is the optimal value 
which ensures that wji(r) approximates well w(r) for 
r < R and simultaneously vanishes rapidly for r > R, 



2 



providing a comparatively fast convergence of the partial- 
wave expansion. The screening radius R must be con- 
siderably larger than the range of the strong interaction 
but from the point of view of scattering theory wr is still 
of short range. Therefore the equations of Ref. [24| be- 
come R dependent. The transition operators WJm where 
a(/3) = 1 and 2 corresponds to initial/final 1 + 3 and 2 + 2 
two-cluster states, respectively, satisfy the symmetrized 
AGS equations 



(G i^G )- 1 P 3 4 - ^34 Ufa G tWGoltfe 



+ U fa Go t ^ Go Ufa , 
Ufa = (G t^Go)- 1 (1-P 34 ) 

+ (l-P 3i )Ufa Got^GoUfa. 
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Here Go is the four free particle Green's function and 
the two-nuclcon t-matrix derived from nuclear potential 
plus screened Coulomb between pp pairs. The operators 



Ufa obtained from 



Ufa — PaG 1 + P a v> R ' Go Ufa, 
Pi = P\2 Piz + -P13 P23, 
P2 = P13 P24, 



(3a) 
(3b) 
(3c) 



are the symmetrized AGS operators for the 1 + (3) and 
(2) + (2) subsystems and Pij is the permutation operator 
of particles i and j. Defining the initial/final 1 + (3) and 
(2) + (2) states with relative two-body momentum p 



l4 fl >(p)> = G i<*>P a |^>(p)>, 



(4) 



the amplitudes for 1 + 3^1 + 3 and 1 + 3^2 + 2 are 
obtained as (p/|T^|pj) = S , a ^(^L H) (p/)|^ ) |^ R) (p i )) 

with Su — 3 and S21 = v3- 

In close analogy with pd elastic scattering, the full scat- 
tering amplitude, when calculated between initial and fi- 
nal p- 3 He states, may be decomposed as follows 



1 (R) — T R 



(R) 



(5) 



with the long-range part t^™ 1 - being the two-body t- 
matrix derived from the screened Coulomb potential of 
the form (TTJ between the proton and the cm. of 3 He, 
and the remaining Coulomb-distorted short-range part 
[Tfa - t% nL ] as demonstrated in Refs. [H Ap- 
plying the renormalization procedure, i.e., multiplying 
both sides of Eq. ([5]) by the renormalization factor ZZ 1 
[lE|20|,inthe-R->oo limit, yields the full 1 + 3 -> 1 + 3 
transition amplitude in the presence of Coulomb 



(P/l^lPi 



(p/l^lPi) 

lim {(p/HT^-trilPi)^ 1 }, 



+ 



(6) 



where the Z R 1 (p f\t° R m - |pj) converges (in general, as a dis- 
tribution) to the exact Coulomb amplitude (p/|ip nL |pi) 
between the proton and the cm. of the 3 He nucleus, 
and therefore is replaced by it. The renormalization fac- 
tor is employed in the partial-wave dependent form as in 

Ref. [H 



Zr 
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with the diverging screened Coulomb p- 3 He phase shift 
rjLR corresponding to standard boundary conditions and 
the proper Coulomb one ax referring to the logarith- 
mically distorted proper Coulomb boundary conditions. 
The second term in Eq. ©, after renormalization by 
Zfi 1 , represents the Coulomb-modified nuclear short- 
range amplitude. It has to be calculated numerically, 
but, due to its short-range nature, the R — > 00 limit is 
reached with sufficient accuracy at finite screening radii 
R. As in pd elastic scattering [15( one needs larger values 
of R for decreasing proton energies, making the conver- 
gence of the results more difficult to reach. Nevertheless 
for E p > 2 MeV the method leads to very precise results 
as we demonstrate in Fig. [T] for the differential cross sec- 
tion da/dfl, proton analyzing power A y , and p- 3 He spin 
correlation coefficient C yy at proton lab energy E p = 4 
MeV. The observables are shown as functions of the cm. 
scattering angle. Fully converged results are obtained 
with R = 12 fm, but already R = 8 and 10 fm results 
are very close to them. The calculations include isospin- 
singlet 27V partial waves with total angular momentum 
I < 4 and isospin-triplet 27V partial waves with orbital 
angular momentum l x < 7, 37V partial waves with spec- 
tator orbital angular momentum l y < 7 and total angular 
momentum J < -t 3 -, 47V partial waves with 1+3 and 2+2 
orbital angular momentum l z < 7 and all initial/final 
p- 3 He states with orbital angular momentum L < 3. The 



charge-dependent (CD) Bonn potential [27J is used. The 
effect of Coulomb is large in the whole angular region, 
particularly for A y where it reduces the magnitude of 
the maximum. The R = fm curve corresponds to the 
so-called Doleschall approximation which clearly fails to 
reproduce the full Coulomb effect. 

In Figs. — H] we compare the results of our calcula- 
tions with data for a number of observables at E„ = 2.25, 
4.0, and 5.54 MeV. In addition to CD Bonn we use AV18 
[2^ , inside- nonlocal outside- Yukawa (INOY04) potential 
by Doleschall 0, [29| , and the one derived from chiral per- 
turbation theory at next-to-next-to-next-to-leading order 
(N3LO) The 3 He binding energy (BE) calculated 

with AV18, N3LO, CD Bonn, and INOY04 potentials is 
6.92, 7.13, 7.26, and 7.73 MeV, respectively; the exper- 



imental value is 7.72 MeV. As in n- 3 H scattering [24(, 



p- 3 He observables depend on the choice of potential; pre- 
dictions with N3LO and AV18 agree best with the cross 
section data but it is INOY04 that provides the high- 
est A y at the peak. If one considers AV18, CD Bonn, 
and INOY04 potentials alone, one might be tempted to 
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FIG. 1: (Color online) Convergence of the p- 3 He scattering 
observables with screening radius R. Results for the differen- 
tial cross section, proton analyzing power Ay, and p- 3 He spin 
correlation coefficient C yy at 4 MeV proton lab energy ob- 
tained with screening radius R = fm (dashed-double-dotted 
curves), 6 fm (dotted curves), 8 fm (dashed-dotted curves), 
f0 fm (double-dashed-dotted curves), and 12 fm (solid curves) 
are compared. Results without Coulomb (dashed curves) are 
given as reference for the size of the Coulomb effect. 



conclude about a possible correlation between observ- 
ables and 3 He BE. Nevertheless, as discussed in Ref. 24), 
N3LO, for reasons not yet fully understood, breaks this 
correlation in the considered energy region. As found in 
Ref. I24| . 4iV S-wave phase shifts correlate with the 3iV 
BE [24J but as the energy increases 4iV P-waves become 



very important as well and behave differently depending 
on the choice of potential. Therefore correlations between 
p- 3 He observables and 3 He BE cannot be established eas- 
ily without further studies, e.g., inclusion of a 3iV force. 

As shown in Figs. [3] — [4] 3 He target analyzing power 
Aq v and p- 3 He spin correlation coefficients Cjk are de- 
scribed quite satisfactorily. This updates the findings of 
Ref. [22| based on AV18 potential where significant dis- 
crepancies were observed for Aq v and C vy . However, the 
proton analyzing power is clearly underestimated by all 
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FIG. 2: (Color online) The differential cross section and pro- 
ton analyzing power A y at 2.25, 4.0, and 5.54 MeV proton 
lab energy. Results including the Coulomb interaction ob- 
tained with potentials CD Bonn (solid curves), AV18 (dashed 
curves), INOY04 (dashed-dotted curves), and N3LO_(dotted 
curves) are compared. The data are from Refs. 0, HD, . 
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FIG. 3: (Color online) 3 He target analyzing power Ao y and 
spin correlation coefficient C yy at 4.0 and 5.54 MeV proton 
lab energy. Curves as in Fig. [5| The data are from Ref. [32I ]. 



potentials. In contrast to low-energy pd elastic scattering 
where variations of the 2N interaction at the maximum 
of A y lead to 10% fluctuations, here we get 15%, which 
means that the 47V system is more sensitive to off-shell 
differences of the 2iV force than the 3N system. 

In Fig. [5] we compare A y for potential INOY04 and 
its version INOY04' @, [5g| with modified 2iV 3 P/ wave 
parameters such that it provides quite satisfactory de- 
scription of A y in low-energy n-d and p-d scattering at 
the cost of being inconsistent with the 2iV data. How- 
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FIG. 4: (Color online) p- 3 He spin correlation coefficients at 
5.54 MeV proton lab energy. Curves as in Fig. [2] The exper- 
imental data are from Ref. [32I ]. 
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FIG. 5: (Color online) Proton analyzing power Ay at 2.25 
and 5.54 MeV proton lab energy. Results for the potentials 
INOY04 (dashed-dotted curves) and INOY04' (solid curves) 
are compared. The experimental data are from Refs. [221, 1 321] . 

ever, for p- 3 He A y disagreement with data still persists. 

Finally, using AV18 potential we investigate the effect 
of 2N magnetic moment (MM) interaction. As for p-d 
scattering [33| it is most visible for A y at low energy 
where at E p = 2.25 MeV it gives rise to a 5.3% increase 
towards the data. At 4 MeV the MM interaction effect 
is reduced to 2.7%. 

In conclusion, we have been able to obtain ab initio 
four-nucleon results for p- 3 He scattering that include the 
Coulomb interaction between the protons for different re- 
alistic local and nonlocal 2N interactions. The reliability 
of the screening and renormalization approach is demon- 
strated. The calculations describe existing data quite 
well except proton A y where there is 25 - 40% discrep- 
ancy at the peak. We find that 47V observables are more 
sensitive than 3N observables to off-shell changes in the 
2N interaction, and that curing A y in low energy 3iV 
scattering through changes in the 2N 3 Pi partial waves 
still gives rise to a p- 3 He A y deficiency. A visible effect 



of 27V magnetic moment interaction is found for A y at 
very low energy. 
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